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We present a generic nonlinear model for current filamentation in semiconductor structures with 
S-shaped current-voltage characteristics. The model accounts for Joule self-heating of a current 
density filament. It is shown that the self- heating leads to a bifurcation from static to traveling 
filament. Filaments start to travel when increase of the lattice temperature has negative impact 
on the cathode-anode transport. Since the impact ionization rate decreases with temperature, this 
occurs for a wide class of semiconductor systems whose bistability is due to the avalanche impact 
ionization. We develop an analytical theory of traveling filaments which reveals the mechanism of 
filament motion, find the condition for bifurcation to traveling filament, and determine the filament 
velocity. 
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I. INTRODUCTION 



Bistable semiconductor systems with S-shaped 
current-voltage characteristics exhibit spontaneous 
formation of current density filaments - high cur- 
rent density domains in the low current density 
environment Current filamentation typically 

develops due to the spatial long wavelength instability 
- known as Ridley instability — of the uniform state 
with negative differential conductance when the device 
is operated via sufficiently large load resistance? 7 ^ 
Fig. 1. Formation of current filaments potentially 
leads to thermal destruction of semiconductor structure 
and thus is an important scenario of semiconductor 
device failure4£ Current filaments have been studied 
in bulk semiconductors)2i£i2ii thin semiconductor 
£i ms ^ii i i2 i i3 4 i4 i i5 anc j i a y erec [ structures of semiconduc- 
tor devices. 16 i 17 i 18 i 19 i 20 i 21 Over last decade, the focus 
has been shifted from static filaments to complex spatio- 
temporal dynamics of current density patternsi 2 ^ 2 ^* 2 ^ 2 ^ 
Modern experimental techniques based on electron 
scanning microscopy^ detection of infrared radiation, 27 
and interferometric mapping 2 — provide means for direct 
observation of current density dynamics during filamen- 
tation. In regard to pattern formation and nonlinear 
phenomena, bistable semiconductors have much in 
common with other spatially distributed active media 
such as gas discharge systems, chemical and biological 
systems, 22 . 29 i 30 i 31 i 32 , 33 

Apart from the well understood bifurcation 
which leads to temporal periodic or chaotic 
oscillationS) 3 ^ 3 ^* 3 ^ 7 . static filament may undergo a 
secondary bifurcation that leads to traveling filament. 
Lateral movement of filaments along the device has 
been observed experimentally for different types of 
semiconductor structures jSLS&SS* 3 ^ Remarkably, in all 



these structures the S-shaped characteristic is associ- 
ated with avalanche impact ionization. Since impact 
ionization coefficients decrease with temperature^ the 
onset of filament motion has been attributed to the 
self-heating of the filament^ - the filament is expected 
to migrate to a colder region as temperature at its 
initial location increases. When current filamentation is 
unavoidable, the migration is desirable in applications 
because the filament motion delocalizes the heating of 
the semiconductor structure and thus reduces the hazard 
of thermal destruction. 

The purpose of this paper is to describe the mechanism 
of filament motion, to develop a generic nonlinear model 
for traveling filament, to find the condition for bifurcation 
to traveling filament, and to determine the filament ve- 
locity. We demonstrate that the transition from static to 
traveling filament due to the Joule self-heating represents 
a generic effect which potentially appears in any bistable 
semiconductor structure provided that transport in the 
cathode-anode direction is sensitive to temperature. 



II. MODEL OF A BISTABLE 
SEMICONDUCTOR SYSTEM 

For many semiconductors and semiconductor de- 
vices current filamentation can be described by a two- 
component reaction-diffusion model which consists of a 
partial differential equation for the bistable element and 
an integro-differential Kirchhoff's equation for the exter- 
nal circuilp 7 *^^ 

Oct 

— = V±{D a (a)V±a) + f(a,u), V_l = e x d x + e y d yi (1) 

d f 

t u — = Uq — u — R I J (a, ujdxdy, t u = RC. (2) 
dt Js 
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Here the variable a(x, y, t) characterizes the internal state 
of the device and the variable u(t) is the voltage over 
the device, Uq is the total applied voltage, J is the cur- 
rent density, R is the load resistance connected in sc- 
ries with the device, C is the effective capacitance of 
the device and the external circuit, and S is the de- 
vice cross-section (Fig. 1). The local kinetic function 
f{a,u) and the J{a,u) dependence contain all informa- 
tion about transport in the cathode-anode (vertical) di- 
rection; the diffusion coefficient D a (a) characterizes lat- 
eral coupling in the spatially extended clement. In the 
bistability range it n < u < u t h (see Fig. 1) the lo- 
cal kinetic function /(a, it) has three zeros a a s,ai n t,a on 
corresponding to off, intermediate, and on branches of 
current voltage characteristic. It has one zero outside 
the bistability range. For the homogeneous state, the 
local dependence a(u) is calculated from /(a, it) = 0, 
and then inserted into J(a, u) in order to determine the 
local current-voltage characteristic J (it) = J(a(u),u). 
Typically d u f,d a j,d u j > 0. The model (JTJ, © be- 
longs to a class of activator-inhibitor models with global 
inhibition^ Variables a and u serve as activator and in- 
hibitor, respectively. The time scale of the activator a 
is r a = <9 Q / -1 ii Specific functional forms /(a,u), j(a, it), 
and D a {a) have been derived for various semiconductors 
and semiconductor structuresi 2 i 7 i 8 i 9 i 10 i 17 i 18 i 20 i 34 i 41 i 42 The 
physical meaning of a depends on the particular type of 
bistable structure: a corresponds to the bias of the emit- 
ter p-n junction for avalanche transistors, 18 ' thyristorsa* 
and thyristorlike structures^ interface charge for het- 
erostructure hot electron diode^i electron charge stored 
in the quantum well for bistable double barrier resonant 
tunneling diode^*^ etc. 

The two-component model Q, is capable to de- 
scribe steady filaments^ as well as temporal oscillation of 
current density filaments i^ 4 *^^ However, its solutions 
do not include traveling filaments. 

The self-heating of the filament has an impact on fil- 
ament dynamics when vertical transport is sensitive to 
temperature. In terms of Eq. £[J it means that the local 
kinetic function / depends on the lattice temperature 
T. We get drf < when temperature suppresses the 
vertical transport. Since the thermal diffusion length £t 
is typically much larger than the device width w, the 
heat dynamics in the device can be modelled by the two- 
dimensional equation 

dT 

cpw — = kw Aj_T + Ju- 7(T - Text), (3) 
at 

Here T corresponds to the mean value of temperature 
along the device vertical direction, c, p, and k are specific 
heat, density, and heat conductivity of the semiconductor 
material, respectively. The second and the third terms on 
the right-hand side of Eq. describe Joule heating and 
cooling due to contact with environment, respectively. 
T ext is the temperature of the external cooling reservoir. 
For multilayer structures with classical transport, T cxt 



is usually a room temperature. The coefficient 7 is a 
heat transfer coefficient (per unit square of the structure) 
that characterizes the efficiency of external cooling. The 
characteristic relaxation time tt, diffusion length £t, and 
propagation velocity vt are given by 

rr^,^ J™, 1*=*= J^£. (4) 
7 V 7 tt y c z p z w 

The modified model is given by the following set of 
equations: 

3d 

— =X7 ± (D a (a)\7 ± a) + f(a,u,T), (5) 
dT 

TT ~dt= 12tA±T + ( Ju/l + Tcxt ~ T ) ' (6) 

du f 
t u — = U -u-RJ J(a,u,T) dxdy. (7) 

The variable T plays a role of the second inhibitor. 

In the following we assume that the device is elongated 
along the x direction (L x 3> L y ) and take only this lateral 
dimension into account. We also neglect direct effect of 
temperature on the current density J in Eq. Q. 

III. CURRENT FILAMENT AND MECHANISM 
OF ITS MOTION 

Current density filament in a long structure represents 
a domain of the high current density state embedded into 
a low current density state (Fig. 2). The width of the 
filament wall If is of the order of £/ ~ \/D a /d a frL For 
T = T e xt the voltage it co over the device with steady 
filament is chosen by the condition known as "equal area 
rule"^ 

/ f(a,u co ,T = T ext )D a (a)da = 0, (8) 

which ensures that production and annihilation of the 
inhibitor a in the filament wall compensate each other. 
When the integral in (JHJ is positive or negative, the bal- 
ance is broken, and the filament walls move in such a 
way that the high current density state or the low cur- 
rent density state, respectively, expand. Since d u f > 0, 
the filament expands for it > it co and shrinks for u < u co . 
The filament has neutral stability with respect to the lat- 
eral shift when T is kept constant £ 

The current-voltage characteristic of the filament in 
long structure is practically vertica&SiSiai and bends 
in the upper and lower parts when the filament be- 
comes narrow (Fig. 1). The bent parts correspond 
to current density intervals [J g; J s + (£ / / L x )J on ] and 
[Jon - {if /L x )J on : J on ], where J on and J g are current 
densities in the high current density state and the low 
current density state, respectively. These intervals are 
negligible for L x > £ f 2ML 

The current filament is stable for sufficiently large load 
resistance R and small relaxation time t u &2& We assume 
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that t u <C T a ,TT, and hence 
sents a constraint imposed on 
by Eqs. 0, ©■ Without loss 
assume that the regime of the 
the current-controlled regime: 
current I sw Uo/Ris constant, 
is directly proportional to the 



Eq. J7J) essentially repre- 
the dynamics determined 
of generality we can also 
external circuit is close to 
Uq u co and the total 
The width of the filament 
total current I 



W = 



II— L X Ly ■ J fi 



off 



J^y won 



where the last equality takes into account that typically 
J on ^> Joft- 

Qualitatively, the mechanism of the filament motions 
in presence of self- heating is the following. With increase 
of temperature the stationary balance JSJ within the fil- 
ament wall is broken due to the temperature dependence 
of the local kinetic function f(a,u,T). As far as the 
temperature profile T(x) is symmetric, the left and right 
walls of the filament are equal and filament would either 
expand or shrink. This is forbidden since the total cur- 
rent is conserved by the global constraint 0). Increase 
of temperature is compensated by deviation of u from 
u co in such a way that the stationary balance within the 
filament wall is restored. Since d u f > and dxf < 0, 
the voltage increases, but the filament stays steady. In 
contrast, for antisymmetric temperature fluctuation the 
balance is disturbed differently in the left and right fila- 
ment walls, becoming positive at one side and negative 
at another side. Potentially, this spontaneous instability 
leads to the motion of a filament as a whole preserving the 
total current. In the traveling filament the temperature 
at the back edge exceeds the temperature at the leading 
edge due to the heat inertia of the semiconductor struc- 
ture. Hence the filament motion becomes self-sustained. 

Below we present an analytical theory of this effect. 
The theory is based on the following assumptions: 

(i) the effect of self-heating is small and can be consid- 
ered as a perturbation, hence the local kinetic function 
can be linearized near T = T ex t and 



f(a,u,T) = f(a,u co ,T = T cxt ) 
+ {u~u co )d u f+{T- 



(10) 



Text) d T /; 



(ii) the transverse dimension of the semiconductor 
structure is large in the sense that L x ^> if,ir, therefore 
only wide filaments with W 3> if are relevant and the 
effect of boundaries is negligible; 

(iii) the width of the filament wall if is much smaller 
than the thermal diffusion length ir and therefore the 
temperature variation within the filament walls can be 
neglected. 



IV. STATIONARY MOTION OF A FILAMENT: 
GENERAL PROPERTIES 

A. Model equations in the co-moving frame 

For stationary motion of a filament with a constant 
velocity v the solution of Eqs. ©, J7J, @ has a form: 

a(x,t) = a(x — vt), T(x,t) — T(x — vt), u — const. (11) 

In the co- moving frame £ = x — vt, Eqs. (0, 0, I© are 

(D a (a)a'y + va' + f(a,u,T) = 0, (12) 
l\ T" + vt t T'+ (Ju/ 7 + T cxt — T) = 0, (13) 
U - u - RL y {J(a,u)) = 0. (14) 

Here the prime (...)' denotes the derivative with respect 
to £ and angular brackets (...) denote integration over £. 
For the sufficiently large L x the boundary conditions are 
given by 



for 



a(0 - a oS (u), T(0 
where 

= Text + Joft (w co )Wco/7 



^±ao, (15) 



(16) 



is the stationary temperature corresponding to the low 
current density state. 



B. Filament velocity and voltage over the structure 

In order to calculate the velocity of the filament and 
the voltage on the structure with traveling filament, we 
note that for given u and T both filament walls can be 
treated as propagating fronts in bistable medium and use 
the standard formula 25 ! 31 1 41 



_ 2j2; { f(a,u 1 T)D a (a)da 
° (D a (a)(a' o y) 



(17) 



where the filament profile is approximated by the sta- 
tionary profile ao(x). Eq. (|17|) implies that the velocity 
v is proportional to the disbalance of production and an- 
nihilation of the activator a in the filament wall. Here 
T is the temperature within the filament wall under con- 
sideration, which is taken as constant according to the 
assumption if <C ir- The factor 2 in the nominator ap- 
pears because ag corresponds to the pattern that consists 
of two fronts. Positive and negative velocities correspond 
to the propagation of the high current density state into 
the low current density state and vice versa, respectively. 

Linearizing the local kinetic function according to Eq. 
((TUj) and taking into account Eq. 10, we obtain 



v{u,T) 



[(« - "co) duf + (T- T C xt) d T f] D a {a) da 



(18) 
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Hereafter, all derivatives are taken at u = u co and T = 

Text- 

Equation (|18|l should be applied separately to the left 
and the right filament walls. In terms of Eq. 1|18|) the 
respective velocities have the same absolute value but 
different signs 



v(u,Tft 



-v(u,T L 



(19) 



where Tr and Tl are temperatures at the right and the 
left walls, respectively. Equations (|18l) and (|19fl yield 
together 



T R - TV 



(D a (a) (a' ) 2 > J ao!f 



d T f D a (a) da, 



,'Tl +Tr 
u = u co + B \ T CKt 



B 



(20) 
(21) 



$a2 dT f Da ^ da 

SaZ d ufD a {a) da 



According to Eqs. (|2(jp. I|21|) the filament velocity is pro- 
portional to the difference of temperatures at the filament 
edges (Tl — Tr), while the voltage deviation from u co is 
proportional to the mean value (T L +T R )/2. For dxf < 
the coefficient B is positive, and hence the voltage in- 
creases with the increase of temperature. The filament 
moves to the right (v > 0) when Tl > Tr and to the left 
(v < 0) when T L < Tr. 

It is worth to mention that formula (|20|l derived for 
the wide filament with narrow walls is a special case of 
the general expression 



(d T fDa(a)T(0a' ) 
(D a (a)(a' o r) 



(22) 



which is applicable for any filament shape and tempera- 
ture profile T(£) satisfying the boundary condition (|15fl . 
Similar to Eq. (|17fl . equation 122|) can be obtained by 
multiplying Eq. i|12|) by D a (a)a' , integrating over £ and 
using the expansion (|10fl . 

Eqs. H2Ufl . (|21|l do not refer exclusively to the case of 
Joule self-heating. They are applicable regardless of the 
origin of the inhomogeneous temperature profile in the 
semiconductor structure. We conclude that when tem- 
perature suppresses the vertical transport (drf < 0), as 
it happens in case of impact ionization mechanism of S- 
type characteristic, the filament generally moves against 
the temperature gradient. In the case when the influence 
of heating is positive (drf > 0), as it happens, for exam- 
ple, when the thermogeneration mechanism is involved, 
the filament moves along the temperature gradient. 



C. Temperature profile in the moving filament 

For a filament with narrow walls (If <C It) the heat 
equation (|13fl is piecewise linear and the temperature 
profile T(£) can be found analytically (see Appendix). 



Consequently, the temperatures Tl and Tr can be pre- 
sented explicitly as functions of the filament velocity v 
and width W. It is convenient to introduce the normal- 
ized difference and the sum of Tl and Tr : 

a _ T L - Tr _ Tl+Tr- 2T+ 

L±LR = T + _ T , ^LR = X* -T ' 

where the temperatures T* and T* are stationary uni- 
form solutions of the heat equation © corresponding to 
uniform on and off states, respectively!*^ 



T* = T cxt + J on (u co ) m co /7, 



(24) 



and T+ is defined by Eq. i|16fl . Typically « T e xt- 
From the solution (|A.2|I of the heat equation we obtain 



A 



LR 



y/l + V' 



Z L R = 1 



(v - cxp (-W\/l+v 2 ^ (25) 
v cosh (Wvj + \J\ + v 2 sinh (W^)]) , 

-j=L= exp (-WVl + v 2 ) (26) 
V 1 + v v ' 



where 



v sinh 



2vt 



(Wvj + y/ 



1 + v 2 cosh 



(Wvj 



W 



W 



The dependencies Alr(v, W) and T,lr{v, W) are shown 
in Fig. 3. The temperature difference is equal to zero for 
v = 0, when the temperature profile is symmetric, and 
reaches maximum at a certain value of v [Fig. 3(a)]. The 
temperature difference decreases with further increase of 
v and vanishes for large velocities (v ^> W/tt) when the 
filament moves too fast to heat the semiconductor struc- 
ture. The average temperature and, according to Eq. 
(|21|l . the voltage u monotonically decrease as v increases 
[Fig. 3(b)]. Hence the onset of filament motion is always 
accompanied by a certain voltage drop. Both (Tl — Tr) 
and (Tl + Tr) increase with filament width W for a given 
v. 



V. SELF-CONSISTENT DETERMINATION OF 
FILAMENT VELOCITY 

In the regime of Joule self-heating, the stationary fil- 
ament motion with a certain velocity occurs when the 
moving filament generates the temperature profile with 
the temperature difference (Tl — Tr) which is exactly 
needed to support this motion. Combining Eqs. I|2(J|I . 
(I23|l . and (|25|l . we obtain a transcendental equation for 
the filament velocity v: 



v 

Vo 
V : - 



A L r(v,W), 

T*-T* 
(D a (a)(a' ) 2 ) 



(27) 

8 T fD a (a)da. (28) 
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Here vq is the upper limit of the filament velocity which is 
achieved for Tl — Tr = T* — T*. This velocity character- 
izes the effect of temperature on the filament dynamics. 
The function Alr is explicitly given by Eq. Al# is 

odd with respect to v, reflecting the symmetry between 
left and right directions of the filament motion. Hence 
nontrivial solutions v* 7^ of Eq. i|27|) always come in 
pairs (v*,—v*). Having this in mind, below we discuss 
only nonnegative solutions v > 0. Note that vq < for 
drf > (positive influence of temperature on vertical 
transport), and therefore Eq. (|28[) has only trivial solu- 
tion v = 0: filament motion is not possible. We shall 
focus on the case drf < when v n > 0. 

It is immediately evident from Fig. 3(a) that Eq. Q27|l 
has either one or two nonnegative roots depending on vq 
and W. Solution v* = 0, corresponding to a steady fila- 
ment, exists for all parameter values. Nontrivial solution 
v* > 0, that corresponds to the traveling filament, exists 
if 



dA 



i'O 



LR 



dv 



> 1. 



(29) 



■u=0 



When the root v = is unique, the corresponding sta- 
tionary filament is stable because Alr(v) < v/vq for 
v > 0. This inequality means that traveling filament gen- 
erates a temperature difference (Tl—Tr) which is not suf- 
ficient to support its motion. The situation changes when 
the condition l|29|l is met: in this case Aijj(u) > v/vo for 
v > 0, and the steady solution is unstable. It means that 
the steady filament looses stability simultaneously with 
appearance of the nontrivial solution v* > which corre- 
sponds to traveling filament. Hence, Eq. H29|l represents 
a condition for spontaneous onset of the filament mo- 
tion. We discuss the bifurcation from static to traveling 
filament in more details in the next section. 

Traveling filament is stable because Alr(v) > v/v 
for v < v* and Alr(v) < v/vq for v > v* . Indeed, these 
inequalities mean that if v decreases the temperature dif- 
ference (Tl — Tr) increases and hence according to i|20|) 
the filament velocity should increase again. In the same 
way, with increase of v the difference (Tl — Tr) decreases 
and therefore the filament slows down. 

In Fig. 4(a) we present numerical solutions v(W) of Eq. 
(|27|l for different values of the parameter ^o. For a given 
vq the motion is possible for filaments whose width ex- 
ceeds a certain threshold Wth(vo) (Fig- 5). With increase 
of W the filament velocity increases and eventually sat- 
urates. We discuss the analytical approximations of the 
front velocity in Sec. VII and Sec. VIII. 

In Fig. 4(b) we present the voltage u in the regime of 
current filamentation obtained by substituting numerical 
solutions of Eq. J2ZJ) into Eqs. The- normalized 

deviation of u from it co 



Am 



,) - (T* - T ext ) 



T*-T+ 



(30) 



is shown, where B is determined by Eq. (|21|) . The voltage 
drop, clearly visible on curves 2-7, is associated with the 



onset of filament motion. With further increase of W the 
voltage increases again. Curve 1 is calculated for vq/vt — 
2.1, which is close to the threshold value vq/vt = 2, and 
practically coincides with the curve for a static filament. 
Since the filament width W is directly proportional to the 
total current / [see Eq. iJSJ], Fig. 4(b) actually represents 
the current-voltage characteristic of the structure with 
traveling filament. 



VI. ONSET OF THE FILAMENT MOTION 

To analyze the onset of filament motion and propaga- 
tion of slow filaments we expand Alr with respect to 
v/vt up to the second order: 




A(W), (31) 



A LR (v,W) 
A(W) 



In this case the temperature profile is close to the sym- 
metric stationary profile given by Eq. (|A.4|) . 

Substituting Eq. (|3*T)l in Eq. (|2*51) , we obtain an explicit 
condition for the onset of filament motion 



(32) 



This condition can be further simplified for narrow and 
wide filaments: 

4^ff) 2>1 fOT W<< ^ ^ 
'" '— + e-3]>l for W~£ T , (34) 



2evT \£t 



^- > 1 for W > £ T , (35) 
2vt 



where e is the natural logarithmic base. Since A(JV) < 1, 
it follows from Eq. I)32|) that regardless to the filament 
width W the filament motion is not possible if v < 2vt- 
If «o > 2vt, the filaments whose width W overcomes 
the threshold Wth determined by Eq. I|32l) start to move, 
whereas smaller filaments remain steady. The depen- 
dence of Wth on Vq/vt is shown in Fig. 5, curve 1. This 
dependence can be approximated as 

W th ^2£ T ,f^ for W^<£t, (36) 

V VQ 

Wth « e T [2e^- + 3 - ej for W th ~£ T , (37) 

Wth »-4r ln ( 1- ^") for W/ th>^T- (38) 

The bifurcation to traveling filament resembles a su- 
percritical pitchfork bifurcation^ at the bifurcation 
point W = Wth the static solution becomes unstable and 
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simultaneously two stable branches corresponding to the 
filaments traveling to the left and to the right appear 
[Fig. 4(a)]. This bifurcation can be understood in terms 
of stability analysis of the current filament performed in 
Ref. |S| for standard two-component model (JIJ , © : The 
spectrum of eigenmodes of a static filament includes a 
neutral mode with zero eigenvalue Xt — which cor- 
responds to translation. Existence of this neutral mode 
reflects the translation invariance of a static filament on 
a large spatial domain. In the extended model JSJ), iJBJ, 
(|7|) the bifurcation to traveling filament is characterized 
by symmetry breaking when Xt becomes positive. This 
becomes possible due to the coupling between the master 
equation JSJ) and the heat equation © when drf < 0. 
Note that for drf > the eigenvalue of the translation 
mode At becomes negative. This corresponds to self- 
pinning of the filament. 



VII. PROPAGATION OF SLOW FILAMENTS 

Substituting Eq. J2H) into Eq. (J23> 

we obtain an ex- 
plicit equation for the filament velocity which is applica- 
ble for slow fronts: 



v(W) = 2v T J—A(W) - 2, 

Vt 



(39) 



where A(W) is defined in Eq. 1|31|) . The asymptotic value 
is given by 



2v T J^ 
vt 



r W 



(40) 



Eq. Ij39|l approximates v(W) with an accuracy of 10% up 
to the parameter value vq/vt = 3 [curve 2 on Fig. 4(a)], 
despite it is derived for v <C «t ■ This wide range of appli- 
cability is due to particulary smooth behavior of Aijj(u) 
to the left of its peak value [see Fig. 3(a)]. 



VIII. PROPAGATION OF FAST FILAMENTS 

In the limit v vt the temperature profile is strongly 
asymmetric [see Eq. (|A.5|) ]. and the temperature at the 
leading edge is close to T*. The expression ((231) can be 
expanded with respect to vt jv : 



A LR (v, W) = 1 - exp 



W 

VTT 



(41) 



/Vr_ 

\ v 



2 — exp 



W 

VTT 



According to Eq. I)A.3|1 the characteristic scale of the tem- 
perature profile for the fast filament is determined by 
(A + ) _1 w vtt which exceeds It- To keep the diffusion 
correction, we expand up to the second order. 

Analytical results are available for the case of a narrow 
filament (W <§; vtt) and a wide filament (W 3> vtt)- 



In the first case (W <C vtt) the temperature inside the 
filament increases linearly, and Tl is much smaller than 
T*. For W/vtt <§C 1, equation (|41|1 reduces to 



W (v T 



(vt 



(42) 



Substituting |g2J| into Eq. O 

we obtain the filament 

velocity 



v(W) 



Wvq ItVt 

2~W 



for 



Tt 



W < vt t . (43) 



Note that in Eq. (|43|) the leading term does not depend 
on the heat diffusion. 

In the case of wide filament (W ^> vtt) the temper- 
ature at the back edge is close the maximum value T*. 
The front velocity is approximated by 



v(W) w v 



2v T 

Vq 



exp 



W \ 

V TT J 

for 



Vq 



2v T 



W > VT T 



(44) 



and saturates at 



2v T 



for 



W 

VqTt 



IX. DISCUSSION 



A. Scales hierarchy - which limiting case is 
relevant? 



The thermal relaxation time tt of the semiconductor 
structure can be straightforwardly determined in exper- 
iment. Therefore it is convenient to use it as basic pa- 
rameter and to express It and vt via tt- Depending 
on the structure width and design, the effective value of 
tt varies from 100-200 ns for transistorlike structures of 
electrostatic discharge protections devices^ (the struc- 
ture width w ~ 10 /lui) to ~ 10-100 ms for power devices 
(the structure width w ~ 100-500 /im)^k It follows from 
Eq. (0} that parameters It and vt are connected to tt 



It — \J Dttt, 



vt 




Dn 



K 



(45) 



where the thermal diffusivity Dt depends only on ma- 
terial parameters. We take Dt — 0.92cm 2 /s and Dt = 
0.25cm 2 /s for Si and GaAs, respectively. Consequently, 
It and vt are of the order of 100 /xm and 10 3 cm/s, re- 
spectively, for small devices (tt ~ 100 ns, w ~ 10 fim). 
We obtain 1 mm and 10 cm/s, respectively, for large 
power devices (jt ~ 10 ms, w ~ 100 /im). Hence in 
most devices the filament width W is of the order of the 
thermal diffusion length It or smaller, and £t is smaller 
but comparable to the transverse dimension of the struc- 
ture L. In particular, it implies that in the regime of 
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self-heating the maximum temperature in the current fil- 
ament is much smaller than T*. For W < it and suffi- 
ciently far from the bifurcation point vq = 2vt the fila- 
ment velocity is well described by the formula (|43[) . where 
the second term can be neglected. 



B. Stimulating the filament motion 

Motion of the filament delocalizes heating of semicon- 
ductor structure and thus is desirable in applications. 
Generally, the start of the filament motion becomes eas- 
ier with increase of vo and decrease of vt and It (see 
Eq. I|32|) and Fig. 5, curve 1). However, these parameters 
cannot be varied independently. Below we focus on the 
effect of the structure parameters which enter our model: 
the heat transfer coefficient 7 and the structure width w. 

Curve 2 in Fig. 5 shows the threshold filament width 
Wth normalized by the quantity (D^tt/vq) 1 ^ 3 which 
does not depend on the heat transfer coefficient 7. Tak- 
ing into account that vt/vq ~ 7 3 ^ 2 , we conclude that de- 
crease of 7 makes the onset of the filament motion easier. 
This occurs due to the increase of vo which, according to 
Eqs. lO, 123}, (HEJ), scales as vq ~ 7 _1 . However, easy 
start of the filament motion in structures with inefficient 
cooling comes at the price of increasing the temperature 
in the filament, which makes static filaments more dan- 
gerous. 

Curve 3 in Fig. 5 shows Wth normalized by the quan- 
tity (Dt/vq) 1 / 2 which does not depend on w. Since 
vt/vq ~ w~ x / 2 , we see that the dependence of Wth on w 
is nonmonotonic. Physically relevant situation W ^ It 
corresponds to the left part of the curve 3 where Wth in- 
creases with w. The effect is due to increase of It ~ w 1 ^ 2 . 
Thus, for the fixed value of vq, thinner structures are 
preferable for the filament motion. 



C. Transient behavior 

Static current filament typically appears due to the 
spatial instability of the uniform state on the middle 
branch of current-voltage characteristic (see Fig. 1). This 
occurs^ on the time scale r a which is typically smaller 
than the thermal time scale tt- First the voltage set- 
tles at u = u co , and only then the temperature starts to 
increase. According to Eq. I|21|) the Joule heating is ac- 
companied by increase of u. As soon as the static filament 
gets hot, the bifurcation to the traveling filament occurs, 
provided that the condition (|29[) is satisfied. The onset of 
filament motion leads to a certain voltage drop, though 
u remains larger than m co [see Eq. (|2"T)l and Fig. 4(b)]. 
The motion can start before the stationary temperature 
profile in the static filament is established, but the non- 
monotonic dynamics of u remains qualitatively the same. 



D. Self-motion and self-pinning 

As it has been pointed out in Sec. VI, we predict, in- 
stead of self-motion, self-pinning of the filament at its 
initial location due to the Joule self-heating in the case 
&rf > 0. For example, it happens when the tempera- 
ture becomes high enough for thermogeneration to set 
in. This typically precedes thermal destruction of the 
semiconductor structure due to local overheating. Eq. 
(|2U[) also suggests an experimental method to distinguish 
between positive (<9t/ > 0) and negative (OtJ < 0) in- 
fluence of temperature on vertical transport by observing 
the filament dynamics in externally applied temperature 
gradient: filaments move along and against the temper- 
ature gradient for 8tJ > and drf < 0, respectively. 

E. Motion of low-current filaments 

S-shaped current- voltage characteristic exhibits formal 
duality between high current density and low current 
density branches. Therefore apart from high-current fil- 
aments, there are patterns in form of low-current fila- 
ments: domains of low current density embedded into 
on state£& Such filaments correspond to the upper part 
of the filamentary current-voltage characteristic (Fig. 1), 
where average current density is close to J on . In the case 
drf < low-current filaments also undergo the bifurca- 
tion from static to traveling filament. Expressions 12011 . 
(|21|l remain valid as they are, whereas in Eqs. (|25|) . i|2ti|) 
T* and T* should be exchanged. The onset of motion re- 
sults in the increase, rather than decrease, of the voltage 
in this case. Duality between on and off states is broken 
when vertical transport is suppressed near the structure 
boundary due to a certain process at the lateral edge of 
the semiconductor structure, e.g., surface recombination 
or surface leak in the p-n junction. (This effect can be 
modelled by Dirichlet boundary conditions a = im- 
posed on the variable a at x = 0, L x £) In this case the 
effect of boundaries can not be neglected for average cur- 
rent densities close to J on even in the limit L x ^£f. The 
current- voltage characteristic, instead of hysteresis, ex- 
hibits a continuous crossover from the filamentary branch 
to the branch of quasiuniform high current density states 
at high current £ Then only high-current filaments are 
observable. 



F. Reaction-diffusion models for traveling spots 

The model (0, ©, J7J belongs to the same class of 
three-component reaction-diffusion models as models for 
trave ling sp ots in active media discussed in Refs. I45H461 
l47u48ll49ll5Ci In contrast to the common two-component 
activator-inhibitor model^^ three-component models 
are capable to describe localized traveling patterns not 
only on a one-dimensional spatial domain, but also on 
spatial domains of higher dimensions^sf The transition 
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from static to traveling spot takes place with increase of 
the relaxation time of the first inhibitor^ 5 - in the same 
way as it occurs in the common two-component model of 
pulse propagation in excitable medial The second fast 
long-range inhibitor plays an essential role only in the 
two-dimensional or three-dimensional case: it prevents 
lateral spreading of the traveling spot which otherwise 
destroys the spatially localized solution and eventually 
leads to development of a spiral wavei^^ 6 . 

This additional inhibition can be either global^ 5 - or 
local4& In the first case the inhibitor has the same value 
in the whole system. This value depends on the mean 
value of other dynamical variables in the system and is 
governed by an integro-differential equation. This cor- 
responds to the global coupling of a spatially extended 
nonlinear system. In the second case the additional in- 
hibitor is governed by a conventional reaction-diffusion 
equation. For traveling spots, the difference between 
these two cases becomes crucial only when several spots 
are considered on a two or higher dimensional domain: 
the system of several spot is unstable when the addi- 
tional inhibition is global, but becomes stable when it 
is local4& This difference vanishes when only one spot 
is present, or the spatial domain is one-dimensional. 46 
Global coupling through the gas phase occurs in the sur- 
face reactions£i*£2a&£i and can be implemented as global 
feedback loop in the light-sensitive Belousov-Zabotinsky 
reactionist. Implementation of the respective control 
loop has allowed to observe localized traveling patterns 
in these systemsi 5 ^ 5 ^ 

In contrast to the three-component models discussed 
in Refs.^lEE in the model ©, ®, both inhibitors 
are needed for the onset of filament motion already in 
the one-dimensional case. Consequently, the bifurcation 
to traveling filament is also different. This reflects the 
fact that we start with a stationary pattern in a bistable 
medium with fast global inhibition (voltage u). This 
global inhibition is due to the external circuit and repre- 
sents an inherent feature of spatio-temporal dynamics of 
a bistable semiconductor: For any evolution of the cur- 
rent density pattern which is accompanied by the varia- 
tion of the total current, the voltage at an external series 
resistance changes, causing the variation of the voltage 
across the device. This inhibition is crucial for the exis- 
tence of current filaments which become unstable when 
the global coupling is eliminated by operating the device 
in the voltage-controlled regime iS&iSi The motion of the 
filament is induced by the effect of another slow diffusive 
inhibitor (temperature T). Similar nonlinear mechanism 
causes the motion of curre nt filam ents obtained by nu- 
merical simulations in Refs. 1231471 However, Refs. 1231471 
focus on a specific type of multilayer thyristorlike struc- 
tures. The model of these devices^ assumes that the 
second inhibitor is a certain internal voltage, not a tem- 
perature. 

It is worth to mention, that self-heating may also trig- 
ger temporal relaxation-type oscillations of a current- 
density filament. This effect has been observed in a re- 



versely biased p -i -n diode and is explained by a similar 
three-component reaction-diffusion models 



X. SUMMARY 

Joule self-heating of a current density filament in a 
bistable semiconductor structure might result in the on- 
set of lateral motion. This occurs when increase of tem- 
perature has negative impact on the vertical (cathode- 
anode) transport. Such negative feedback takes place 
when bistability of semiconductor structure is related 
to the avalanche impact ionization, because the im- 
pact ionization rate decreases with temperature^ Ex- 
amples of such devices are avalanche transistors^ re- 
versely biased p -i -n diodes in the regime of avalanche 
injection, electrostatic discharge protection devices oper- 
ated in the avalanche regime^ multilayer thyristorlike 
structuresi2ii21 Traveling filaments can be consistently 
described by a generic nonlinear model ©, ©, ©■ 

Generally filaments move against the temperature gra- 
dient with a velocity proportional to the temperature dif- 
ference at the filament edges [Eq. 1)20(1] . In the regime of 
Joule self-heating the strength of coupling between the 
master equation JSJl, which controls the current density 
dynamics, and the heat equation © can be characterized 
by single parameter vq [Eq. (|28|l ]. which has a dimen- 
sion of velocity. The filament velocity v is determined 
by the transcendental equation (|27|) . The condition for 
the spontaneous onset of filament motion l|29|) depends 
on the ratio of vq and the thermal velocity vt, and on 
the ratio of the filament width W and the thermal dif- 
fusion length if. Filament motion is never possible for 
vq < 2vt- For vq > 2vt, static filaments whose width 
W exceeds a certain threshold Wth [Eq. I|32ll . curve 1 in 
Fig. 5] become unstable and start to travel. The filament 
velocity v and the voltage on the structure u are shown 
in Fig. 4. For most semiconductor structures W < It, 
and sufficiently far from the bifurcation point vq = 2vt 
the filament velocity can be approximated by a truncated 
version of Eq. I|43|l 

Heating of the static filament is accompanied by an in- 
crease of the voltage u [Eq. (|2ip]. With onset of the 
filament motion the voltage drops, but remains higher 
than the voltage u co [Fig. 4(b)]. 

Our analytical theory does not cover narrow filaments, 
when the flat top of current density profile does not exist 
or the filament width W is comparable to the width of 
the filament wall if , as well as cylindrical filaments which 
appear when lateral dimensions of semiconductor struc- 
ture L x and L y are comparable. However, the bifurcation 
from static to traveling filament and the mechanism of 
filament motion remain qualitatively the same in these 
cases. 
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\2v T J 2v T 



Eqs. (|A.2|I allow to determine the temperature difference 
(T L + T R ) and the mean temperature (T L + T R )/2 [Eqs. 
(|23|1 . I|25|l . (|26|) ] which are needed for self-consistent cal- 
culation of the front velocity v and the voltage u. 



APPENDIX: SOLUTIONS OF THE HEAT 
EQUATION 

Assuming that the front walls are thin, we present Eq. 
(G3J) as (see Fig. 2): 



t\ T" + vt t T' + {T* -T) = 



(A.l) 



W , W 

for -y <?< T' 



l\T" + vt t T' + (T*-T) = 

W W 
for £< — -and£>y. 

Here the middle of the filament is at £ = 0. T* and T* are 
defined by Eqs. I|24l) . The solution of this piecewise 
linear equation is given by 



T(0 = T*-(T*-T* 
x 



(A.2) 



sinh ( — r— I exp(A £) 



A+ - X- 



for 



_ W 



no = t* - g f: ( a ( .xp 



A+ - A" 



A- £ 



W 

Y 



A exp 



A+ £ 



T 



for —<£<—. 
2 s 2 

T(£) = T*-(T*-T+) 

2A" . , /A+HA „ , 

x aT-a^^H^J ^ 



for 



£< 



where A + and A are eigenvalues of Eq. (|A.1|) : 



Steady temperature profile. — For v = we get A + = 
— A~ = 1/£t an d the corresponding symmetric tempera- 
ture profile is given by 



for £ > 



IT 



(A.4) 



T ( £) = T*- ( T*-T.) exp cosh (|- 

W ^ FT 

for T <e< T' 

T ( £) = T* + (T*-T.) sinh (|-) exp (A 
IT 

for £ < -y- 



Temperature profile in the fast filament. — For n » we 
get A + = (vtt) , A - = — oo. The corresponding profile 
is strongly asymmetric: 



W 



T(£) = T* for £> 

T(£) = T*-(T*-T*)exp 

W W 
for -— < £ < — , 
2 s 2 



(A.5) 



2£- PF 
2vtt 



T* + 2(T* - T*) sinh f ) exp f -1- 

\ 2vtt J \ vtt 



for £ < 



if 



i 



\ 2vt / 2ut 



(A.3) 



Electronic mail: rodin@physik.tu-berlin.de 

1 B. K. Ridley, Proc. Phys. Soc. London 82, 954 (1963). 

2 A. F. Volkov and Sh. M. Kogan, Sov. Phys. Usp. 11, 881 
(1969) [Usp. Phys. Nauk 96, 633 (1968)]. 

3 V. L. Bonch-Bruevich, I. P. Zvyagin, and A. G. Mironov, 



Domain Electrical Instabilities in Semiconductors (Consul- 
tant Bureau, New York, 1975). 

4 S. M. Sze, Physics of Semiconductor Devices (Wiley, New 
York, 1981). 

5 E. Scholl, N onequilibrium Phase Transitions in Semicon- 



10 



ductors (Springer- Verlag, Berlin Heidelberg, 1987). 

6 M. Shaw, V. Mitin, E. Scholl, and H. Grubin, The Physics 
of Instabilites in Solid State Electron Devices (Plenum 
Press, New York, 1992). 

7 A. Wacker and E. Scholl, J. Appl. Phys. 78, 7352 (1995). 

8 A. Alekseev, S. Bose, P. Rodin, and E. Scholl, Phys. Rev. E 
57, 2640 (1998). 

9 A. F. Volkov and Sh. M. Kogan, Sov. Phys. JETP 25, 1095 
(1967) [Zh. Eksp. Teor. Fiz. 52, 1647 (1967)]. 

10 F.G. Bass, V. S. Bochkov, Yu. Gurevich, Sov. Phys. JETP 
31, 972 (1970) [Zh. Eksp. Teor. Phys. 58, 1814 (1970)]. 

11 V. Novak, C. Wimmer, and W. Prettl, Phys. Rev. B 52, 
9023 (1995). 

12 F.-J. Niedernostheide, J. Hirschinger, W. Prettl, V. Novak, 
and H. Kostial, Phys. Rev. B 58, 4454 (1998). 

13 V. Novak, J. Hirschinger, F.-J. Niedernostheide, W. Prettl, 
M. Cukr, and J. Oswald, Phys. Rev. B 58, 13099 (1998). 

14 J. Hirschinger, F.-J. Niedernostheide, W. Prettl, and V. 
Novak, Phys. Rev. B 61, 1952 (2000). 

15 G. Schwarz, C. Lehmann, and E. Scholl, Phys. Rev. B 61, 
10194 (2000). 

16 I. V. Varlamov, V. V. Osipov, and E. A. Poltoratsky, Sov. 
Phys. Semicond. 3, 978 (1970) [Fiz. Tekh. Poluprovodn. 3, 
1162 (1969)]. 

17 V. V. Osipov and V.A. Kholodnov, Sov. Phys. Semicond. 
4, 1033 (1971) [Fiz. Techn. Poluprovodn. 4, 1216 (1970)]. 

18 V. V. Osipov and V.A. Kholodnov, Mikroelectronika 2, 
529 (1973) (In Russian). 

19 D. Jager, H. Baumann, and R. Symanczyk, Phys. Lett. A 
117, 141 (1986). 

2(1 A. V. Gorbatyuk, P. B. Rodin, Solid-State Electron. 35, 
1359 (1992). 

21 A. V. Gorbatyuk and F.-J. Niedernostheide, Phys. Rev. B 
59, 13157 (1999); 65, 245318 (2002). 

22 B. S. Kerner, V. V. Osipov, Autosolitons (Kluwer Aca- 
demic Publishers, Dordrecht, 1994). 

23 Nonlinear Dynamics and Pattern Formation in Semi- 
conductors and Devices, edited by F.-J. Niedernostheide 
(Spinger- Verlag, Berlin, 1995). 

24 K. Aoki, Nonlinear Dynamics and Chaos of Semiconduc- 
tors (Institute of Physics Publishing, Bristol, 2000). 

25 E. Scholl, Nonlinear Spatio- Temporal Dynamics and Chaos 
in Semiconductors (Cambridge University Press, Cam- 
bridge, 2001). 

26 A. Wierschem, F.-J. Niedernostheide, A. Gorbatyuk, and 
H.-G. Purwins, Scanning 17, 106 (1995). 

27 F.-J. Niedernostheide, B. S. Kerner, and H.-G. Purwins, 
Phys. Rev. B 46, 7559 (1992). 

28 D. Pogany, S. Bychikhin, M. Litzenberger, E. Gornik, G. 
Groos, and M. Stecher, Applied Physics Letters 81, 2881 
(2002). 

29 M. S. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

Y. Kuramoto, Chemical Oscillation, Waves and Turbu- 
lence (Springer- Verlag, Berlin, 1988). 
31 A. S. Mikhailov, Foundations of Synergetics (Springer- 
Verlag, Berlin, 1994). 



32 Self-organization in Activator- Inhibitor Systems: Semi- 
conductors, Gas Discharge, and Chemical Active Media, 
edited by H. Engel, F.-J. Niedernostheide, H.-G. Purwins, 
E. Scholl (Wissenshaft&Technik, Berlin, 1996). 

33 Evolution of Spontaneous Structures in Dissipative Con- 
tinious Systems, edited by F. H. Busse and S. C. Miiller 
(Springer, Berlin, 1998). 

34 A. Wacker and E. Scholl, Z. Phys. B: Condens. Matter 93, 
431 (1994). 

35 F.-J. Niedernostheide, H. Schulze, S. Bose, A. Wacker, and 
E. Scholl, Phys. Rev. E 54, 1253 (1996). 

36 S. Bose, P. Rodin, and E. Scholl, Phys. Rev. E 62, 1778 
(2000). 

37 F. Plenge, P. Rodin, E. Scholl, and K. Krischer, 
Phys. Rev. E 64, 056229 (2001). 

38 K. Penner, Journal de Physique: Colloque C4 49(9), 797 
(1988). 

39 D. Pogany, S. Bychikhin, E. Gornik, M. Denison, N. 
Jensen, G. Groos, and M. Stecher, Proc. IEEE Int. Re- 
liab. Phys. Symp., IRPS2003, Dallas, Texas, 2003, p.241. 

40 G. Wachutka, IEEE Transactions on Electron Devices 38, 
1516 (1991). 

41 M. Meixner, P. Rodin, E. Scholl and A. Wacker, Eur. Phys. 
J. B 13, 157 (2000). 

42 P. Rodin and E. Scholl, J. Appl. Phys. 93, 6347 (2003). 

43 Here we negelect dependence of T* and T* on u, assuming 
u — u co . This does not change the results. 

44 J. Guckenheimer and P. Holmes, Nonlinear Oscillations, 
Dynamical Systems, and Bifurcations of Vector Fields, 
Vol. 42 of Applied Mathematical Sciences (Springer- Verlag, 
Berlin, 1983). 

45 K. Krischer and A. Mikhailov, Phys. Rev. Let. 73, 3165 
(1994). 

46 CP. Schenk, M. Or-Guil, M. Bode and H.-G. Purwins, 
Phys. Rev. Let. 78, 3781 (1997). 

47 F.-J. Niedernostheide, M. Or-Guil, M. Kleinkes, and H.-G. 
Purwins, Phys. Rev. E 55, 4107 (1997). 

48 H. Hempel, I. Schebesh, and L. Schimansky-Geier, Eur. 
Phys. J. 2, 399 (1998). 

49 M. Or-Guil, M. Bode, CP. Schenk, and H.-G. Purwins, 
Phys. Rev. E 57, 6432 (1998). 

50 M. Bode, A.W. Liehr, CP. Schenk, and H.-G. Purwins, 
Physica D 161, 45 (2002). 

51 F. Mertens, R. Imbihl, and A. Mikhailov, J. Chem. Phys. 
101, 9903 (1994). 

52 M. Falcke and H. Engel, J. Chem. Phys. 101, 6255 (1994). 

53 M. Bertram and A. Mikhailov, Phys. Rev. E 67, 036207 
(2003). 

54 M. Bertram, C. Beta, M. Pollmann, A. Mikhailov, H. 
Potermund, and G. Ertl, Phys. Rev. E 67, 036208 (2003). 

55 A.N. Zaikin and A.M. Zabotinsky, Nature 225, 535 (1970). 

56 E. Mihaliuk, T. Sakurai, F. Chirila, and K. Showalter, 
Phys. Rev. E 65, 065602(R) (2002). 

57 B. Datsko, Semiconductors 31, 146 (1997) [Fiz. Tekh. 
Poluprovodn. 11, 250 (1996)]. 



11 




FIG. 1: Current-voltage characteristic of a bistable structure. The average current density I/S is shown, where S = L x L y 
is the cross-section of the structure. Unstable middle branch with negative differential conductance is depicted by the dashed 
line. The hold and threshold voltages are denoted as «h and w t h, respectively. The vertical branch at u — u co corresponds to 
a static filament. The inset shows a sketch of a bistable semiconductor structure operated in the external circuit with load 
resistance R, capacitance C, and bias Uo- 




FIG. 2: Current density profile in a filament (thick line). Thin line denotes the temperature profile in the static filament. 
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FIG. 4: The filament velocity v (a) and normalized voltage on the structure Au (b) as functions of the filament width W 
for different values of the parameter vq. Am is defined by Eq. (J3UJ. Curves 1 to 7 correspond to v = 2.1, 3, 10, 20, 50, 100, 200, 
respectively. Peaks on the voltage curves are related to onset of the filament motion. 
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FIG. 5: Normalized threshold filament width Wth corresponding to onset of the filament motion as a function of vt/vo- 
Curve 1 shows Wth normalized to the thermal diffusion length: Wo = It- Curve 2 shows Wth normalized to the quantity 
Wo = (D^tt/vo) 1 ^ 3 , which does not depend on the heat transfer coefficient 7. Curve 3 shows W t h normalized to the quantity 
Wo = \J Dt I vq , which does not depend on the structure width w. 



